Astronomy <fc Astrophysics manuscript no. 2724 


February 2, 2008 


(DOI: will be inserted by hand later) 





Interstellar deuterated ammonia: From NH 3 to ND 3 

E. Roueff 1 , D.C. Lis 2 , F.F.S. van der Tak 3 , M. Gerin 4 and P.F. Goldsmith 5 



in 
o 
o 

(N 
Sh 

Oh 
< 
O 



> 

in 

o 
in 
o 

Oh 
6 



13 



Laboratoire Univers et Theorie, UMR 8102 du CNRS, Observatoire de Paris, Section de Meudon, Place Jules 
Janssen, 92195 Meudon, France 

2 Downs Laboratory of Physics 320-47, California Institute of Technology, Pasadena, CA 91125, USA 

3 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 

4 Laboratoire d'Etude du Rayonnement et de la Matiere en Astrophysique, UMR 8112 du CNRS, Observatoire 
de Paris and Ecole Normale Superieure, 24 Rue Lhomond, 75231 Paris Cedex 05, France 

5 Department of Astronomy and National Astronomy and Ionosphere Center, Cornell University, Ithaca, NY 
14853, USA 



Received / Accepted 

Abstract. We use spectra and maps of NH2D, ND2H, and ND3, obtained with the CSO, IRAM 30m and Arecibo 
telescopes, to study deuteration processes in dense cores. The data include the first detection of the hyperfine 
structure in ND2H. The emission of NH2D and ND3 does not seem to peak at the positions of the embedded 
protostars, but instead at offset positions, where outflow interactions may occur. A constant ammonia fractionation 
ratio in star-forming regions is generally assumed to be consistent with an origin on dust grains. However, in the 
pre-stellar cores studied here, the fractionation varies significantly when going from NH3 to ND3. We present a 
steady state model of the gas-phase chemistry for these sources, which includes passive depletion onto dust grains 
and multiply saturated deuterated species up to five deuterium atoms (e.g. CD^~). The observed column density 
ratios of all four ammonia isotopologues are reproduced within a factor of 3 for a gas temperature of 10 K. We 
also predict that deuterium fractionation remains significant at temperatures up to about 20 K. ND and NHD, 
which have rotational transitions in the submillimeter domain are predicted to be abundant. 

Key words. ISM: molecules - Molecular processes - Stars: circumstellar matter - Stars: formation 



1. Introduction 

In recent years, the chemistry of dense cores has 
been revised with the almost simultaneous discov- 



ery of multiply deuterated molecules (D?CO: iTurnerl 


Il99l ICeccarelli et al. 


1998: ND,H: iRoueff et all l200fi 


ND^:lLis et al.ll2002bl 


van der Tak et alJl20()3 CHD,OH. 


CD^OH: IParise et al. 


l2002l 120041) on the one hand. 



letions of CO on the other hand 
fact, the two processes 



and of very large dep 

ijBacmann et alJ l2002|) . In tact, tne two processes are 
closely related since the condensation of CO and other 
abundant gas phase molecules favors deuterium fraction- 
ation, by decr easing the destructio n rate of deuterated 
molecular ions ijBacmann et alJ2003j) . As multiply deuter- 
ated molecules are found in dense cores prior to star for- 
mation (as in Barnard 1 and LDN 1689N), or in very 
young protostars (as in NGC 1333 IRAS 4A), they can 
be used for studying and identifying the earliest stages of 
star formation. While an analysis of the line profiles gives 
clues about the gas dynamics in the envelope surrounding 
the protostar, or in the dense core, the molecular abun- 
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dance and especially the fractionation provides informa- 
tion on the physical conditions (temperature, depletion) 
and on the chemistry. It is therefore of interest to improve 
our understanding of deuterium chemistry. For this pur- 
pose, accurate fractionation ratios are needed for multiply 
deuterated molecules. 

Two main paths are invoked for understanding the 
large deuterium fractionation observed in dense cores. The 
first is based on gas phase chemistry and invokes the ion- 
molecule deuterium exchange reactio ns taking plac e at 
low temperatures, as first proposed bv lWatsonf l|l974h. In 
this framework, deuterated molecules are built in succes- 
sive reactions starting with eithe r H2D 4 " (and also D2H+ 
and D+), CH,D+ or CHD+ llRoberts fc Millarl l200fl 
iRoberts et ail 120031 iGerlich et alJ I2002I) . The confirmed 
detections of H 2 D+ in the envelopes of disks encompassing 
the protostars NGC 1333 IRAS 4A and IRAS 16293-2422 
l|Stark et al.lll999UStark et aljl2004h and in the pre-stellar 
core LDN 1544 IjCaselli et al.ll2003l) . as well as the recen t 
detection of D 2 H+ in LDN 1689N l)Vastel et al.ll2004t) . 
have provided strong observational support for this theory. 
The second path for forming deuterated molecules is based 
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on grain chemistry l)Tielenslll983|) . As hydrogenation on 
grain mantles has been proposed t o explain the large 
abundances of solid phase CH^ OH llDartois et al.1 [l999, 
iGibb et al"ll2000l IPontoppidan et alJl2003[) . the analogous 
process with deuterium would build deuterated molecules 
on grains. Because in dense cores, the gas phase abundance 
ratio of atomic D/H is probably orders of magnitude larger 
than the elemental abundance ratio of 1.5 x 10 -5 , the sur- 
face chemistry leads also to large deuterium fractionation, 
which scales as the (gas-phase) D/H ratio 1 . While the sur- 
face processes seem interesting, there is no secure confir- 
mation of the role of dust mantles for deuterium frac- 
tionation, as neither HDO nor deute rated methanol has 
yet been detected in the solid phase ijDartois et al.ll20o"3 
IParise et al]l2003j) . 

These two paths may lead to different frac- 
tionation ratios for mu ltiply deuterated molecules 
i Rodgers fc Charnlevl l200lh . because the formation of 
deuterated molecules involves largely different processes 
and chemical reactions. For gas phase processes, the frac- 
tionation ratio at each stage results from competing chem- 
ical reactions (fractionation reactions substituting deu- 
terium atoms for hydrogen atoms, destruction reactions 
with other molecules and dissociative recombination re- 
actions removing the fractionation, etc.). The fractiona- 
tion ratio must be calculated with a comprehensive chem- 
ical model and is expected to be different for the different 
stages of deuterium fractionation of a given species, since 
the specific reactions arc different. In a simple model of 
grain surface chemistry, the abundance ratios of various 
deuterated species are expected to be similar, since each 
such ratio is proportional to the gas phase D/H abun- 
dance ratio. For all bonds with a heavy atom, the proba- 
bility of bonding with a H or D atom scales with the deu- 
terium/hydrogen ratio. All bonds behave independently 
in this process. 

Four different isotopologues of ammonia, including 
fully deuterated ammonia, ND3, as well as NH2D and 
ND2H, can be observed in dense cores. The ammonia 
data can thus provide the information required for under- 
standing the formation of multiply deuterated molecules 
in dense cores. This paper presents new data on deuter- 
ated ammonia in a small sample of dense cores, some of 
which are forming stars. The observations are presented 
in Sect. |3 and the results described in Sect. Sect. 0] 
presents comprehensive gas phase chemical models aimed 
at understanding the roles of gas density, depletion, and 
temperature in the fractionation of ammonia. We summa- 
rize our results in Sect.0 

2. Observations 

The rotational transitions and frequencies of the vari- 
ous isotopologues of ammonia that have been observed 
are shown in Table ^ The frequencies refer to the 

1 This paper uses full names (e.g., deuterium) to denote el- 
ements, and symbols (e.g., D) to denote atomic species. 



strongest hyperfine compon ent. For Einstein A— values 
and statistical we ights, see iTine et all ll200fil) (NH^D) 
iRoueff et aljEffltl fNDoffl a,ndlvan der Tak et rtlbmA 
andlLis et a.lJ l|20n2rl (ND 3 ). 



Table 1. Parameters of observed rotational transitions 



Species 


Transition 


E u 


V 








K 


MHz 


NH 2 D 


In 


-loi 


20.7 


85926.3 


NH 2 D 


In 


-loi 


21.3 


110153.6 


ND 2 H 


ho 


- loi 


18.7 


110812.9 


ND 2 H 


lio 


-loi 


18.4 


110896.7 


ND 2 H 


2 2 o 


-2ia 


56.8 


206971.9 


NH 2 D 


3 22 


— 3i 2 


120. 


216562.6 


ND 3 


loi 


-Ooo 


14.9 


309909.7 



2.1. Millimeter-wave observations 

Observations of NH2D and ND2H lines near 110, and 
220 GHz were carried out at the 30-m telescope of the 
Institut de Radio Astronomie Millimctrique (IRAM) on 
Pico Veleta, Spain, in August 2002. We also use addi- 
tional 86 GHz observations performed in November 2001 
and September 2002. The front ends were the facility re- 
ceivers A100, B100, A230 and B230, and the back end was 
the Versatile Spectral and Polarimetric Array (VESPA) 
autocorrelator. During night time, the weather was good, 
with «4 mm of precipitable water vapour and T sys Ril40 K 
at 110 GHz, while during daytime, the water column and 
T sys roughly doubled. At 110 GHz, the telescope has a 
FWHM beam size of 22". The pointing of the telescope 
was stable to within 2 arcsec. Data were calibrated onto a 
T m b scale by multiplying by 1.25, the ratio of the forward 
coupling efficiency to the mai n beam efficiency. 

Sources were selected from lShah fc Woottenl l)200l|) to 
have high N H?D column densiti es, but the positions used 
here are f roml.Iiiina et alJ (1999) for Orion and Monoceros 
and from lMcMidlin et al.l 1 20001 ) for Serpens. Table[2]lists 
source positions and LSR velocities. Only upper limits 
were obtained toward NGC 2264G and CB 17. 

To locate the 'deuterium peak' in each source, maps 
were made of the NH2D In — > loi emission near 110 GHz. 
These maps include up to 18 points, at half-beam (11") 
spacing. In the case of NGC 1333, a 55-point map was 
made, encompassing the locations of both previous ND3 
detections. The integration time was 2 minutes per point, 
using frequency switching with a frequency shift of 8 MHz. 
Subsequently, deep integrations on the ND2H lines were 
performed at the locations of the NH2D peaks. 

2.2. Submillimeter-wave observations 

Submillimeter observations were carried out in 2002 and 
2003 using the 345 GHz facility receiver and spectrome- 
ters of the Caltech Submillimeter Observatory (CSO) on 
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Table 2. Positions and velocities of the observed sources. 



Source 


R.A. 


Dec. 


Vlsr 




J2000 


J2000 


km s _1 


NGC 1333 


03:29:10.3 


+31:13:32.2 


7.0 


Barnard 1 


03:33:20.8 


+31:07:34.8 


6.8 


LDN 1544C 


05:04:16.6 


+25:10:47.8 


7.4 


OMC 2 


05:35:26.3 


-05:09:49.4 


11.2 


HH 1 


05:36:17.3 


-06:46:08.2 


9.2 


LDN 1630 


05:46:07.4 


-00:13:41.8 


10.0 


NGC 2264C 


06:41:11.6 


+09:29:25.1 


8.0 


LDN 134N 


15:54:08.5 


-02:52:01.0 


2.4 


LDN 1689N 


16:32:29.5 


-24:28:52.6 


3.8 


LDN 483 


18:17:29.7 


-04:39:38.3 


5.5 


S 68N 


18:29:47.5 


+01:16:51.4 


8.5 


NGC 2264G 


06:41:11.0 


+09:56:00.8 


4.6 


CB 17 


04:04:38.0 


+56:56:11.0 


-4.6 



Mauna Kea (Hawaii). The data were taken under aver- 
age weather conditions (225 GHz zenith opacity ~0.1). 
The CSO FWHM beam size at 309 GHz is -25", and 
the main beam efficiency at the time of the observations, 
determined from total-power observations of planets, was 
~60%. Typical calibration uncertainties are ^25%. The 
pointing of the telescope was stable to within 5 arcsec. 
We used the 1.5 GHz, 500 MHz, and 50 MHz bandwidth 
acousto-optical facility spectrometers. The spectral reso- 
lution of the high-resolution, 1024 channel 50 MHz spec- 
trometer is ~0.15 MHz or 0.15 kms -1 . All spectra were 
taken in the double beam switched observing mode with 
a secondary chopper throw of 180". 



The main beam plus first sidelobe efficiency of the Arecibo 
Telescope is ~70% at 1.6 GHz. 

3. Results 

3.1. Maps of NH2D emission 



20 



- • 



-20 - • 



20 -20 

R.A- offset (arcsec) 

Fig. 1. Map of f T m bdV of the main hyperfine compo- 
nent of the NH2D In ^ lot (ortho) 85.9 GHz transition 
in Barnard 1. The coordinates of the (0,0) position are 
R.A.(J2000) = 03:33:20.8, Dec.(J2000) = +31:07:34.8 . 
Greyscale steps are: 0.5 to 5.0 by 0.5 K km s _1 . 



2.3. Centimeter-wave observations 

Observations of the ND3 inversion lines in the Barnard 1 
cloud were carried out between 2003 October 31 and 
November 11 at the Arecibo Observatory in Puerto Rico, 
using the new L-Band Wide facility receiver. Four sub- 
correlators of the interim spectrometer system were used, 
each having 1024 lags. A high-resolution subcorrelator 
with a total bandwidth of 0.781 MHz (0.144 km s" 1 
velocity resolution) was centered on the frequency of 
the strongest ND3 hyperfine component, 1589.0178 MHz 
(Tabled) ■ Two additional subcorrelators with bandwidths 
of 3.125 MHz, centered on the frequencies of 1588.4043 
and 1588.6316 MHz, covered the remaining ND3 hyperfine 
components with a velocity resolution of 0.576 km s _1 . In 
addition, a 1.563 MHz bandwidth subcorrelator, centered 
on the frequency of the 1667.3590 MHz OH line, was used 
to monitor system performance during the long integra- 
tion. Quasars B02231+313 and B0340+048 were used for 
pointing and calibration. The data were taken in position 
switching mode with the reference position at a right as- 
cension offset by 5 minutes of time. The typical system 
temperature referred to above the earth's atmosphere was 
~25 K. The total on source integration time was 756 min. 



Figures t° El show the spatial distribution of NH2D 
emission as observed with the IRAM 30m. In Barnard 1 
and HH-1, the emission has a distinct peak at a slight 
offset from our nominal central position given in Table 
|2J On the other hand, the emission in OMC-2, S68N and 
NGC 2264C is not strongly peaked. The map of NGC 1333 
indicates that NH2D does not peak at the position of the 
protostar IRAS 4A at (0,0). Instead, the emission peaks in 
two lobes, one located 10-20" N of IRAS 4A and the other 
to the W-SW of it, roughly between offsets (20,0) and 
(30,-20). This distribution follows that of the DCO + 3^2 
emission, mapped at the CSO bv lLis et alJ l(2004h . and is 
consistent with the observation that the ND3 emission in 
NGC 1333 is stronger at the offset p osition (23,-6) than 
at IRAS 4 A l|van der Tak et al.ll2002tl . 

Table 01 reports the column density of NH2D at those 
positions where the fitting of the hyperfine s t ructur e 
was successful, calculated following iTine et all l|200f](L 
Emission of NH2D is clearly detected in other cases, as dis- 
played in the maps, but no optical depth could be derived 
as the individual spectra have a low signal to noise ratio. 
The HFS method within the CLASS package gives the to- 
tal optical depth of the line and the radiation temperature, 
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R.A. offset (arcsec) 



R.A. offset (arcsec) 



Fig. 2. Map of the integrated NH 2 D In — l i (para) 
110.2 GHz emission in HH-1. The coordinates of the 
(0,0) position are R.A.(J2000) = 05:36:17.3, Dec.(J2000) 
= -06:46:08.2 . Greyscale steps are 0.3 to 0.6 by 0.15 
K km s _1 . Unfilled circles denote positions where no emis- 
sion was detected. 




R.A. offset (arcsec) 



Fig. 3. As previous figure, for OMC-2. The coordi- 
nates of the (0,0) position are R.A.(J2000) = 05:35:26.3, 
Dec.(J2000) = -05:09:49.4 . Greyscale steps are 0.1 to 0.5 
by 0.1 K km s -1 . 



assuming a beam filling factor of unity. Since this latter as- 
sumption may not always be fulfilled, the table also gives 
column density estimates for the cases T ex =5 K and 10 K. 
We calculated the partition functions at these two con- 
sidered temperatures. The partition functions of NH 2 D, 
ND2H, and ND3 for an excitation temperature of 10 K 
are 9.37, 9.51, and 38.3, respectively. The corresponding 
values for T cx —5 K are 4.47, 4.21, and 17.23, respectively. 

The absence of an NH 2 D peak at NGC 1333 IRAS 4A 
could partly be due to excitation effects. However, the 
3 22 — * 3i 2 line of NH 2 D is not detected at this position. 
The upper limits for this line in NGC 1333, Barnard 1 and 
NGC 2264C and G are T mb <70-100 mK in 0.11 km s" 1 



Fig. 4. As previous figure, for S68N. The coordi- 
nates of the (0,0) position are R.A.(J2000) = 18:29:47.5, 
Dec.(J2000) = +01:16:51.4 . Greyscale steps are 0.1 to 0.7 
by 0.1 K km s -1 . 




R.A. offset (arcsec) 



Fig. 5. As previous figure, for NGC 2264C. The coordi- 
nates of the (0,0) position are R.A.(J2000) = 06:41:11.6, 
Dec.(J2000) = +09:29:25.1 . Greyscale steps are 0.1 to 0.7 
by 0.1 K km s -1 . 



channels. For S68N, OMC-2, HH-1 and CB 17, which were 
observed during daytime, the corresponding upper limits 
are T mb <200-300 mK. 

3.2. Spectra of ND2H emission 

Emission in the lio — > lot lines of ND 2 H was detected 
only in Barnard 1; upper limits were obtained in other 
sources. The spectra have rms levels of T m b=5-ll mK in 
0.11 km s" 1 channels. The ND 2 H lines toward Barnard 1 
are the strongest observed to date, and the spectra (Fig.[7J) 
are the first ones to show the hyperfine components of the 
transition. The hyperfine structure has been calculated by 
L. Coudert (priv. comm.) and the corresponding frequen- 
cies and relative strengths are given in Table 0] Our fit 
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Table 3. Observed NH2D line parameters and derived column densities. Numbers in parentheses denote uncertainties 
in units of the last decimal. 



oource 


OllocL 


1 - 
> 


AT/ 


TV 


T 


T 

± cx 




iv 






£LrCS6C 


-1 

Km s 


-1 

Km s 


K 
l\ 




K 

IV 




If) 14 nm~ 

iu cm 


-2 


Barnard 1: 




















+0 


+0 


6.78(1) 


0.76(2) 


1.3(2) 


4.0(5) 


4(3) 


5.8(7) 


4.3(5) 


2.9(3) 


+11 


+11 


6.76(1) 


0.69(4) 


2.1(5) 


2.5(5) 


5(3) 


2.1(5) 


2.4(5) 


1.6(3) 


+11 


+0 


6.82(1) 


0.67(2) 


2.2(2) 


2.7(3) 


6(3) 


2.2(2) 


2.6(2) 


1.7(2) 


+11 


-11 


6.88(1) 


0.72(3) 


1.9(4) 


2.6(5) 


5(3) 


2.6(5) 


2.7(5) 


1.8(3) 


+0 


-11 


6.80(1) 


0.70(3) 


1.6(4) 


2.7(6) 


5(3) 


2.9(6) 


2.7(5) 


1.8(4) 


-11 


-11 


6.66(3) 


0.97(6) 


1.1(4) 


2.5(8) 


4(3) 


5(2) 


3(1) 


2.3(8) 


-11 


+11 


6.78(2) 


0.74(5) 


1.1(3) 


3.4(8) 


4(3) 


5(1) 


3.6(9) 


2.4(6) 


+22 


+0 


6.83(2) 


0.73(6) 


2(1) 


1.2(7) 


6(5) 


0.9(5) 


1.3(7) 


0.9(5) 


+22 


+11 


6.72(2) 


0.61(5) 


2(1) 


1.6(7) 


6(4) 


1.0(4) 


1.4(6) 


0.9(4) 


+22 


-11 


6.90(2) 


0.70(6) 


1.2(5) 


2.4(9) 


4(3) 


3(1) 


2.4(9) 


1.6(6) 


+0 


-22 


9.21(4) 


0.47(7) 


0.8(6) 


3(2) 


4(4) 


4(3) 


2(1) 


1(1) 



a Values in column 8 are for T ex equal to the estimate in column 7, and in columns 9 and 10 for T CJ 
Uncertainties in N include only those arising from the uncertainties in r. 



=5 K and 10 K, respectively. 




Table 4. Frequencies of the hyperfine components of the 
lio — * lot transitions of para and ortho ND2H . 



R.A. offset (arcsec) 



Fig. 6. As previous figure, for NGC 1333. The coordi- 
nates of the (0,0) position are R.A.(J2000) = 03:29:10.3, 
Dec.(J2000) = +31:13:32.2 . Greyscale steps are 0.1 to 0.9 
by 0.1 K km s" 1 . 



to the hyperfine structure, indicates a low optical depth 
(r < 0.1) for the lines. 

The column densities of ND2H have been calculated 
from the o bserved line strength s, using spectroscopic con- 
stants from lRoueff et all I 2000l h We assume optically thin 
emission, but take deviations from the Rayleigh- Jeans law 
into account. Table presents the results. We assume T cx 
= 10 K for positions close to young stellar objects and 5 K 
for other positions. The num bers given in the table sup er- 
sede previous values given in Ivan der Tak et al.l ([2003) . 

Simultaneous observations of the ND2H 220 — 2n line 
did not result in any detections. Upper limits to T m b range 



Transition 


z;(MHz) 


Relative Strength 


Ortho 


F=l-1 


110811.1130 


1/12 


F=l-2 


110811.7265 


5/36 


F=2-l 


110812.2972 


5/36 


F=l-0 


110812.6467 


1/9 


F=2-2 


110812.9107 


5/12 


F=0-1 


110814.0735 


1/9 


Para 


F=l-1 


110895.0298 


1/12 


F=l-2 


110895.6433 


5/36 


F=2-l 


110896.2108 


5/36 


F=l-0 


110896.5635 


1/9 


F=2-2 


110896.8243 


5/12 


F=0-1 


110897.9823 


1/9 



from 11 to 22 mK in 0.11 km s _1 channels (T sys of 400- 



500 K), except for S68N and NGC 1333 (33,-22), where 
worse weather (T sys ~1000 K) limited the sensitivity to 
40-50 mK. 

3.3. Spectra and maps of A/D3 emission 

Table H3 and Figure |S1 summarize the results of the ND3 
observations. In addition, Figure El shows a map of ND3 
in LDN 1689N, with overlaid contours of DCO+ (3-2) 
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Table 5. Observed ND2H line parameters (or la upper limits) and derived column densities. Numbers in parentheses 
denote uncertainties in units of the last decimal. Positions are the central positions given in Table [21 unless an offset 
is indicated. 



• • 

Position 


ort no / 




A 1/ 

L\ V 


J JmbU' 


T 

J - cx 


1\ (o+p ) 




para 


km s 


km s 


niK km s 


IV 


1U UI1 


Barnard 1 





6.76(2) 


0.67(5) 


64(4) 


5.0 


97(6) 




P 


6.97(2) 


0.34(5) 
0.5 b 


30(4) 


5.0 


75(10) 


HH-1 







< 2.7 


10.0 


< 1-1 




p 




0.5 b 


< 2.2 


10.0 


< 1.5 


S68N 







0.8 6 


< 2.8 


10.0 


< 1.1 




p 




0.8 6 


< 2.8 


10.0 


< 1.9 


NGC 1333: 














(0,0) 







1.6° 


< 3.6 


10.0 


< 1.5 




P 




1.6 a 


< 3.6 


10.0 


< 2.5 


(0,11) 







1.6 a 


< 2.2 


10.0 


< 0.9 




P 




1.6 a 


< 2.2 


10.0 


< 1.5 


(23,-6) 







1.0 a 


< 2.1 


5.0 


< 3.2 




P 




1.0 a 


< 1.8 


5.0 


< 4.5 


(33,-22) 







1.0 a 


< 3.5 


5.0 


< 5.3 




P 




1.0 a 


< 3.9 


5.0 


< 9.7 



a Value taken from ND 3 Ivan der Tak et al-lEoO^ 
b Value taken from NH 2 D (Table 

Table 6. Observations of ND3. Offset positions from nominal central position are given in arc seconds, enclosed in 
parentheses. 



Position 


Vo 


AV 


/ T mh dV 


Tqx. 


N 




km s" 1 


km s _1 


K km s" 1 


K 


10 12 cm" 2 


NGC 1333 


7.2 


1.6(6) 


0.071(21) 


10 


0.29 (9) a 


NGC 1333 (23,-6) 


6.5 


1.0(2) 


0.144(43) 


5 


1.3(4) 


NGC 1333 (-24,106) 


7.8 


0.7(2) 


0.078(17) 


5 


0.70(15) 
2.8(7) b 


Barnard 1 (0,0) 


6.5 


0.82(7) 


0.307(19) 


5 


Barnard 1 (-31,113) 






< 0.045 


5 


< 0.4 


LDN 1544C 


6.9 


0.34(10) 


0.120(24) 


5 


1.1(2) 


LDN 1630 (20,35) 






< 0.096 


5 


< 0.9 


LDN 134N 






< 0.036 


5 


< 0.32 


LDN 1689N (0,0) 


3.3 


0.42(4) 


0.165(17) 


5 


1.49(15) 


LDN 1689N (-10,-10) 


3.4 


0.44(3) 


0.202(23) 


5 


1.8(2) 


LDN 1689N (0,-20) 


3.3 


0.42(3) 


0.288(26) 


5 


2.6(2) 


LDN 483 


5.3 


1.5(3) 


0.081(19) 


5 


0.7(2) 



a van der Tak et al. (2002) ; b Lis et al. (2002). The column density given here is given for T cx = 5 K. 



and 1.3 mm dust continuum emission. The emission is 
seen not to peak on IRAS 16293, ~ 90" to the east of 
the map center, but rather at the outflow interaction 
region between the YSO and a pre-stellar core. Many 
deuterated molecules peak at this position, possibly be- 
cause the outflows from IRAS 16293 A and B, which are 
not co-eval, have desorbed grain mantles and compressed 
the gas, which h as subsequently cooled l|Lis et alJl2002al 
IStark et alJ20ol . Indeed, the emission peaks at a velocity 
offset from that of the core (Table 0). 



3.4. A/D3 inversion lines 

The high-resolution spectrum at the frequency of the 
strongest hyperfme component of the 1.59 GHz ND3 inver- 
sion line in Barnard 1 is shown in Figure lTUI No emission is 
detected with a 3cr upper limit for the integrated intensity 
of 0.011 Kkm s _1 , computed over a 2 km s^ 1 wide velocity 
interval and corrected for the main beam efficiency. This 
upper limit has to be further corrected for the source cou- 
pling to the 2'.5 Arecibo beam. The spatial distribution 
of ND3 in the Barnard 1 cloud has not been determined. 
Therefore to estimate the coupling efficiency, we assume 
that the ND 3 emission has the same spatial distribution 
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Fig. 7. Spectra of the ND 2 H 1 10 - lot ortho (bottom) 
and para (top) transitions towards Barnard 1, with our fit 
to the hyperfine structure overplotted. 



as the DCO + ( 3-2) emission, as mapped with the CSO 
l|Lis et al.11200^) . The DCO+ intensity in a 2.5' beam is 
decreased by a factor of ^2 compared to the peak value 
in the 30" CSO beam. We thus estimate a 3tx upper limit 
for the integrated intensity of the 1589.0178 MHz ND 3 
inversion line in the CSO beam to be ^0.022 Kkm s _1 . 

To compute the corresponding upper lim i t for th e ND3 
column density we use Eq. (2) of iLis et al.l (|2002bfl . with 
A = 5.23 x 10~ n s"\ g u = 24, and E u = 11.97 K, The 
ND3 partition function Q = 38.3 for T ex — 10 K, leading 
to a 3c upper limit of 3.0 x 10 13 cm~ 2 for the ND 3 col- 
umn density, after including a factor of 2 correction for 
the remaining hyperfine components (see Table 0). For 
T ex — 5 K, the partition function Q = 17.2 and the cor- 
responding 3cr upper limit for the ND3 column density is 
7.0 x 10 13 cm -2 . These li mits are an improvement over 
previous Effelsberg results ijvan der Tak et af1l2003[) . but 
still factors of ~20 above the v alues derived from the ro- 
tational lines {Lis et, a.l]l2002hh . 



Table 7. ND3 in version frequencies within t he Jk = li 
level derived from van Veldhoven etTaH (|2002|) . 



Transition 


i/(MHz) 


Relative Strength 


F=0-1 


1587.4645 


1/9 


F=2-l 


1588.3969 


5/36 


F=l-1 


1589.0114 


1/12 


F=2-2 


1589.0086 


5/12 


F=l-2 


1589.6231 


5/36 


F=l-0 


1590.5545 


1/9 




(0,0) 



(0,-20) 



-1 1 1 J 1 I 1 

LDN1689N n (-10,-10) 



f- 

LDN134N (x3) 
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Fig. 8. Spectra of ND3 309.9 GHz emission, observed with 
the CSO. The Barnard 1 (Bl) spectrum is from Lis et al. 
(2002) and the NGC 1333 IRAS 4A spectrum from van 
der Tak et al. (2002). 



3.5. Summary of ammonia deuterium fractionation 

Table El lists the col umn densities and fractionation ratios 
for ammonia isotopologues in our sample. Only sources 
where ND 3 , or three other isotopologues have been de- 
tected are listed. Care was taken that the data are in 
matching beams (20-40") whenever possible. For example, 
there are signific ant discrepancies in the colum n density 
deriva tions from Sha h fc Woottenl l|200l|) and Hatchcll 
l|2003h towards NGC 1333 IRAS 4A. The NH 3 data for 
NGC 1333, Barnard 1 and LDN 134N have been taken 
with the Effelsberg telescope with a ~40" beam. Data for 
LDN 1544C and LDN 483 refer to a larger beam ( ~100") 
(see references listed in Table |5J), and the NH3 column 
density for LDN 1544 and LDN 483 may thus be under- 
estimated. 

Ammonia fractionation of the various isotopologues is 
seen to vary from source to source, whereas grain sur- 
face chemistry predicts that the deuterium fractionation 
should scale as atomic D/H ratio, and be constant to 
within a factor of a few. However, we cannot rule out 



8 



Roueff et al: From NH 3 to ND 3 . 



20 



3 -20 



-40 




0.3 



-DCO + (3-2) 1.3mm cont 

J i i i I i i i I i 



0.2 



0.1 



20 



-20 



Aa (arcsec) 



Fig. 9. Intensity of the ND 3 emission in LDN 1689N inte- 
grated over velocities between 2.8 and 5 km s _1 , (greyscale 
and white contours), with overlayed contours of DCO + (3- 
2) emission (solid black lines) and 1.3 mm dust contin- 
uum emission (dashed black lines). The coordinates of the 
(0,0) position are R.A.(J2000) = 16:32:29.5, Dec.(J2000) 
= -24:28:52.6 . Contour levels are 50, 70, and 90% of the 
peak for ND3, 50 to 90% of the peak, with an interval 
of 10% for DCO+, and 4, 6, and 8% of the peak (toward 
IRAS 16293, outside the plot) for dust continuum. 




Fig. 10. High-resolution spectrum at the frequency of the 
1589.0178 MHz ND3 hyperfine component in Barnard 1 
(left scale). The dotted line shows the 1667.3590 MHz OH 
line at the same position (right scale). 

the grain s urface scenario, as the pr edictions are only 
qualitative ijRodeers fc Charnlevl200l|) . For most sources, 
[NH 2 D]/[NH 3 ] and [ND 2 H]/[NH 2 D] ratios are of the same 
order of magnitude, ~ 10 - 20%, while the [ND 3 ]/[ND 2 H] 
ratio is significantly smaller (~ 2 — 5%). The assumption 
of one single excitation temperature for all lines in one 



source may not be valid, since the lines have very differ- 
ent excitation conditions (Table ^J. This uncertainty may 
affect the column density ratios in Table |S1 because the 
column densities of some isotopologues are derived from 
ground state lines and others from excited states. From 
the observational side, the largest uncertainty is thus in 
the radiative transfer modelling of the data. Collisional 
cross sections for ammonia isotopologues are needed for a 
more accurate analysis. 

The ND 3 to NH 3 ratio varies between - 1 - 2 x 10~ 3 
in LDN 1689N, Barnard 1, and the depleted core of 
LDN 1544, where H ?D + and D?H+ have been detected 
jCaselli et al.ll2003l and IVastel et alJl2004[) . to less than 
1.6x 10~ 4 for the dark cloud LDN 134N. Such values repre- 
sent prodigious enhancement factors over the cosmic ratio 
of 12 orders of magnitude to less than 11 orders of magni- 
tude which deserve specific chemical studies. We present 
in the next section a gas phase chemical model and the 
corresponding detailed results. 

4. A model of gas-phase deuterium chemistry 

4.1. Approach and general results 

We have updated the preliminary gas phase chemical net- 
work, which was used to interpret ou r first detection o f 
ND 3 towards Barnard 1, as presented in lLis et"al] 1 2002fJl . 
We have introduced the formation of H 2 , HD, and D 2 
on grain surfaces by f ollowing the simple approach pre- 
sented in iLe Bourlot et ail l)l995[) . The role of the grains 
(passive, although very significant) is to trap the condens- 
able species, which are then no longer available in the gas 
phase. The resulting chemical abundances may differ con- 
siderably from the so-called typical values, relevant for a 
standard molecular cloud such as TMC-1. Such effects are 
labelled as depletion, and observations show a clear cor- 
relation of deuterium fractionation and CO depletion in 
dense pre-stellar cores. This trend results, to first order, 
from the shorter collision time between dust grains and 
molecules when the density increases. 

In order to mimic this effect, we consider three different 
models, in which the density and depletion are varied. We 
consider a "standard" molecular cloud similar to TMC1 
with a density of n(H 2 ) = 10 4 cm -3 , a temperature T of 
10 K and abundances relative to molecular hydrogen of 
TMC-1, i.e. [C]/[H 2 ] = 7x 10~ 5 ; [0]/[H 2 ] = 2 x 10~ 4 ; 
[N]/[H 2 ] = 2 x 10~ 5 . The relative abundance of sulfur is 
not well known and is taken as [S]/[H 2 ] = 3.6 x 10~ 7 . In 
addition, we introduce a representative metal [M] which 
does not contribute to the molecular complexity, and only 
undergoes charge transfer and radiative recombination re- 
actions, with [M]/[H 2 ] = 3 x 10~ 8 (low metal case). This 
model is referred as model 1. For all models, we neglect 
possible deuterium depletion on grains and have adopted 
a gas phase deuterium abundance of |D]/[H 1 =1.5 x 10~ 5 
as derived from the latest determinations llLinskvl 12003). 
Models 2 and 3 have densities n(H 2 ) = 10 5 cm" 3 and 10 6 
cm~ 3 and carbon and oxygen depletion factors of 5 and 15, 
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JUUI Lt 


IN 113 


NH2D 


ND2H 








[INU3J / 






10 15 cm" 2 


10 14 cm" 2 


10 13 cm" 2 


10 12 cm" 2 


[NH 3 ] 


[NH 2 D] 


[ND 2 H] 


[NH3] 


NGC 1333 


1.4±0.4 a 


3.9±0.8 a 


<0.5 h 


0.29±0.09 d 


0.28±0.14 


<0.01 


>0.06 


2.1±lxl0" 4 


Barnard 1 


2.5±0.2 ,; 


5.8±0.7 9 


8.6±0.8 h 


2.8±0.7 / 


0.23±0.05 


0.15±0.03 


0.033±0.01 


1.1±0.5 xlO -3 


LDN 1544C 


0.20±0.035 fc 


0.26±0.01 fe 




1.1 ±0.2 e 


0.13±0.03 






5.5±lxl0~ 3 


LDN 134N 


2.0±0.2 J 


2.0±0.4 6 


1.0±0.2 C 


<0.32 e 


0.10±0.03 


0.05±0.02 


< 0.032 


<1.6xl0~ 4 


LDN 1689N 


1.8±0.2 fc 


3.4±0.7 ; 


7.6±1.1 ; 


1.8±0.2 e * 


0.19±0.05 


0.22±0.05 


0.024±0.01 


1.0±0.5 xl0~ 3 


LDN 483 


1.4±0.1 m 






0.7±0.2 e 








5.0±1.2xl0" 4 



a lHatchelll fa00$) . b iTine et alJ ]2000|L c iRoueff et all J2000h . 20 % accuracy assumed, d Ivan der Tak et all <2002lL e Table El 
* the - 10" , -10" offset value has been given as it corresp onds to the ND2H peak obs erved bv iLoin ar^Tt ^lTlpOulTv . f iLis et all 
l|2002bh . 9 Table El h Table 151 ' | Bachiller etafl l|l990li . 10 % accuracy assumed, j llJngerechts et all (jl98(Tl . 10 % accuracy 
assumed, fc IShah fe Woottenl feOOllLl O % acc uracy assumed, ' lLoinard et alJ <l200ll) . reanalysis of the para transition of NH2D 
from the original data. m Ijiiina et al] lll999l) . lAnglada et al] ]l997h . 10% accuracy assumed. 



respectively. Nitrogen is kept constant in the three models 
as ISO spectra have not indicated that ammonia ices are 
abundant. The temperature has been kept equal to 10 K 
and the cosmic-ray ionization rate is taken as 2 x 10~ 17 
s _1 . Table lists the predicted molecular abundances for 
the three models. 

We consider 210 gas phase species, including H 2 , H 2 0, 
CH 4 , NH 3 , H 2 S, H 2 CS, all deuterated variants thereof, 
and corresponding ions. A total of 3000 chemical reac- 
tions are included in the chemical network. In addition 
to the reactions of Hj, CH;j", C 2 H 2 f with HD, w hich are 
drivin g the deuteration as already introduced in IWatsonl 
we have considered the subsequent deuteration of 
molecular hydrogen by introducing the reactions of H 2 D + 
and D 2 H + with HD, D, a nd D? as studied experimen - 
tally at low temperatures bv lGerlich fc Schlemmerl i|2002lb 
Some of these react ions hav e also been introduced by 
iRoberts et all l|2004h and by IWalmslev et~aTl l|2004l) for 
the completely depicted case. Table ITU1 gives the fraction- 
ation reactions, the rate coefficients of the forward cor- 
responding reaction and the endothermicities involved in 
the reverse reactions. These values are obtained from the 
ground level energies without differentiating between or- 
tho and para forms of the various species involved. 

As already discussed in ILis et all l|2002b|) and iRouefll 
(200^), we have taken special care regarding the dissocia- 
tive recombination reactions and introduced several up- 
dates such as the new dissociative r ecombination rate 
coefficient of H 3 + l|McCall et aljl2003jb the new branch- 



Table 10. Fractionation reactions involving H3 , 
and C 2 H+ with HD, D 2 and D. 



CH^ 



ing ratios of t h e diss ociative recombination of N 2 H + by 
iGeppert et alJ ^2004) , who found that NH is a major 



product channel an d the recent measurement of NDj 
llOiekull et all 120041) . In addition, we have generalized 
the result of recent studies of dissociative recombination 
reactions of partially deuterated molecular ions (H 2 D + , 
HDO + and HD 2 + ) which show that the release of H 
atoms is favoured compared to that of deuterium atoms. 
We assume that the branching ratios differ by a factor of 2 
between H and D ejection, following the recommendations 





k 


AE 


Ref. 




cm 3 s _1 


K 


Ref. 


H+ + HD # H 2 D+ + H 2 


3.5(-10) 


232 


1 


H 2 D+ + HD ^ D 2 H+ + H 2 


2.6(-10) 


187 


1 


H+ + D 2 ^ H 2 D+ + HD 


1.75(-10) 


152 


1 


H+ + D 2 ^ D 2 H+ + H 2 


1.75(-10) 


340 


1 


D 2 H+ + HD # D+ + H 2 


2.00(-10) 


234 


1 


H 2 D+ + D 2 ^ D 2 H+ + HD 


3.50(-10) 


109 


1 


D 2 H+ + D 2 ^ D+ + HD 


3.50(-10) 


100 


2 


H+ + D ^ H 2 D+ + H 


1.0(-9) 


632 


1 


H 2 D+ + D # D 2 H+ + H 


1.0(-9) 


600 


1 


D 2 H+ + D # D+ + H 


1.0(-9) 


600 


4 


CH+ + HD # CH 2 D+ + H 2 


2.6(-10) 


375 


1 


CH+ + D 2 # CD 2 H+ + H 2 


4.4(-10) 


375 


5 


CH+ + D 2 # CH 2 D+ + HD 


6.6(-10) 


375 


1 


CH 2 D+ + HD ^ CD 2 H+ + H 2 


1.0(-9) 


370 


1 


CH 2 D+ + D 2 ^ CD 3 + + H 2 


3.0(-10) 


370 


3 


CH 2 D+ + D 2 ^ CD 2 H+ + HD 


9.0(-10) 


370 


3 


CD 2 H+ + HD ^ CD 3 + + H 2 


5.64(-10) 


370 


3 


CD 2 H+ + D 2 ^ CD 3 + + HD 


4.5(-10) 


370 


3 


C 2 H+ + HD # C 2 HD+ + H 2 


7.5(-10) 


550 


1 


C 2 HD+ + HD # C 2 D+ + H 2 


7.5(-10) 


300 


1 


C 2 H+ + D 2 ^ C 2 HD+ + HD 


7.0(-10) 


550 


2 


C 2 HD+ + D 2 ^ C 2 D+ + HD 


7.0(-10) 


550 


2 


C 2 H+ + D # C 2 HD+ + H 


7.5(-10) 


250 


2 


C 2 HD+ + D ^ C 2 D+ + H 


7.5(-10) 


150 


2 



s: l.lGerlich et al 


J2002I) 


Walmsle^t^l ( 


2004); 



); 2. estima t e: 3. | Smith et all 
5. lAnicichl J2003I) 



of M. Jensen (private communication) 2 . The influence of 
this assumption is discussed in Section fOl 

We can derive the following general trends from these 
models: 

1. The fractional ionization decreases as the density in- 
creases, as expected. The main source of electrons is 



In lLis et al] (l2002bf) we had assumed a factor 3 to 1 
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Table 9. Prediction of fractional abundances relative to molecular hydrogen 



Model 1 Model 2 Model 3 

n(H 2 ) (cirr 3 ) 10 4 10 s 10 6 



Depletion 1 5 15 



x( e -r 


6.05 


X 


10" 


-8 


3.74 


X 


10" 


-8 


2.91 


X 


10" 


-8 


H 


2.25 


X 


10" 


-4 


2.26 


X 


10" 


-5 


2.27 


X 


10" 


-b 


HD 


2.75 


X 


10" 


-5 


2.65 


X 


10" 


-5 


2.56 


X 


10" 


-5 


D 


1.22 


X 


10" 


-6 


4.43 


X 


10" 


-7 


8.12 


X 


10" 


-8 


D 2 


5.84 


X 


10" 


-7 


1.49 


X 


10" 


-6 


2.18 


X 


10" 


-6 


CO 


6.79 


X 


10" 


-5 


1.33 


X 


10" 


-5 


4.47 


X 


10" 


-6 


c 


1.72 


X 


10" 


-6 


5.64 


X 


10" 


-7 


1.64 


X 


10" 


-7 


H 3 + 


1.03 


X 


10" 


-8 


3.01 


X 


10" 


-9 


4.78 


X 


10" 


-10 


H 2 D+ 


5.80 


X 


10" 


-10 


4.45 


X 


10" 


-10 


1.05 


X 


10" 


-10 


D 2 H+ 


3.13 


X 


10" 


-11 


6.47 


X 


10" 


-11 


2.33 


X 


10" 


-11 


D 3 + 


1.29 


X 


10" 


-12 


8.35 


X 


10" 


-12 


5.20 


X 


10" 


-12 


NH 


4.89 


X 


10" 


-y 


6.80 


X 


10" 


-y 


1.66 


X 


10" 


-y 


ND 


1.18 


X 


10" 


-9 


4.67 


X 


10" 


-9 


1.70 


X 


10" 


-9 


NH 2 


6.01 


X 


10" 


-8 


8.48 


X 


10" 


-8 


1.25 


X 


10" 


-8 


NHD 


2.99 


X 


10" 


-9 


7.64 


X 


10" 


-9 


1.51 


X 


10" 


-9 


ND 2 


3.05 


X 


10" 


-11 


2.05 


X 


10" 


-10 


5.09 


X 


10" 


-11 


NH 3 


5.49 


X 


10" 


-8 


8.82 


X 


10" 


-8 


2.00 


X 


10" 


-8 


NH 2 D 


2.71 


X 


10" 


-9 


8.80 


X 


10" 


-9 


2.94 


X 


10" 


-9 


ND 2 H 


8.82 


X 


10" 


-11 


6.10 


X 


10" 


-10 


2.37 


X 


10" 


-10 


ND 3 


3.44 


X 


10" 


-12 


5.14 


X 


10" 


-11 


2.14 


X 


10" 


-11 


HCN 


7.13 


X 


10" 


-8 


3.78 


X 


10" 


-8 


7.92 


X 


10" 


-9 


DCN 


2.40 


X 


10" 


-9 


3.15 


X 


10" 


-9 


8.90 


X 


10" 


-10 


HNC 


3.29 


X 


10" 


-8 


1.95 


X 


10" 


-8 


4.02 


X 


10" 


-9 


DNC 


1.54 


X 


10" 


-9 


2.35 


X 


10" 


-9 


6.28 


X 


10" 


-10 


HCO+ 


9.67 


X 


10" 


-9 


1.02 


X 


10" 


-9 


8.01 


X 


10" 


-11 


DCO+ 


2.85 


X 


10" 


-10 


6.85 


X 


10" 


-11 


8.15 


X 


10" 


-12 


N 2 H+ 


2.66 


X 


10" 


-10 


5.28 


X 


10" 


-10 


1.68 


X 


10" 


-10 


N 2 D+ 


9.24 


X 


10" 


-12 


4.25 


X 


10" 


-11 


1.80 


X 


10" 


-11 



a fractional ionization. 



provided by the representative metal which may only 
recombine or charge transfer and does not contribute 
to the chemical complexity. 

2. The deuterium reservoir is provided by HD. However, 
as the density increases, molecular D2 becomes a signif- 
icant deuterium containing species. It is then impor- 
tant to account properly for the gas phase reactions 
involving doubly deuterated hydrogen. 

3. Model 2 presents the highest molecular fractional 
abundances for the N-bearing species, while HCN, 
HNC, and HCO + follow the same pattern as C and 
CO and have higher fractional abundances in model 1. 

4. The gas phase chemistry changes significantly between 
models 2 and 3. Indeed, with the high depletion as- 
sumed for model 3, the Hjj" chemistry is driven by re- 
actions with HD rather than reactions with the abun- 
dant neutrals O, CO, etc... This results in significant 
variations in the fractional abundances of saturated 
species such as NH 3 , HCN, HNC, and HCO+ and 
N2H" 1 " molecular ions. 

5. D2 is the sole species, which has higher abundance in 
model 3, as compared to model 2. 



6. The combination of high density and depletion results 
in very efficient production of deuterated species. The 
abundances of the multiply deuterated species such as 
D+ and ND 3 , are enhanced by between 10 and more 
than 12 orders of magnitude relative to those expected 
from the cosmic deuterium abundance, in the various 
models we have considered. 

7. All models predict high fractional abundances of am- 
monia progenitors such as NH, ND, NH2, NHD and 
ND 2 . 

4.2. Comparison with the observations 

We now focus on the comparison of the deuterium frac- 
tionation values of the various species with the differ- 
ent observed values. Table ITD reports the deuterium frac- 
tionation ratio obtained with the three previously de- 
scribed models and the observed values towards a sam- 
ple of clouds. We also report the fractionation ratios of 
nitrogen containing molecules calculated when statistical 
approximation is assumed for the reaction products com- 
ing from dissociative recombination reactions, when the 
latter have not been studied in the laboratory. The result- 
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Model 2 


Model 3 
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DCN : HCN 
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0.019 


0.049 
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DNC : HNC* 
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DCO+ : HCO+ 
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0.18 d 


0.08 e 


0.04 h 


N 2 D+ : N 2 H+ 


0.035 
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0.11 


0.35 d 
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* models employing a statistical approximation for the reaction products in dissociative recombination reactions of non studied 

d euterated mo l ecular ions m ^ ^ 

11 IVastel et alJ ll2004h . para-D 2 H + / ortho-H 2 D + b Tabled c iTurneil (1200 ll) based on N(HN 13 C) x60, d lTine et al] ll200Ch. e 
iLis et , al] l|2002ah . based on N(H 13 CN) and N(H 13 CO+) x60. f iHirota et al] l|200ll) . based on NfHN 13 C) x60. 9 Iffirota et al] 
j20na|L based on N(HN 13 C) x60, h ICaselli et all (|2002h . i lOerin et al] (|200lTl 



ing fractionation ratios differ at most by a factor of 2 in ated formaldehyde bv ICeccarelli et all l|2002f) in the same 
the nitrogen containing species. pre-stellar core. 



Following previous remarks, the deuterium fractiona- 
tion ratio is seen to increase from model 1 to model 3. The 
D 2 H+ / H 2 D+ ratio found in LDN 1689N involves signif- 
icantly excited levels (the upper level energies of H2D + 
and D2H 4 " involved are 104 and 83 K, respectively) and 
may not reflect the actual val ue of the full D9H 4 " / H ^D + 
ratio calculated in o ur models. IWalmslev et al] l)2004h and 
iFlower et alJ ( 20A4tl have considered a completely depleted 
chemistry with specific para and ortho species and find a 
p-D2H + / o-H2D + consistent with the observations. 

However, the ensemble of data for this pre-stellar core 
(sometimes also called IRAS 16293E) is spectacular. The 
level of deuteration is particularly high. lStark et alJ l)2004|) 
deduce a depletion of CO of the order of 5 compared to 
TMC-1 (from [C 18 0] / [H 2 ]= 3 x 1CT 8 and [CO] : [C 18 0] 
= 500 : 1), equal to the value assumed in model 2, and 
a temperature in the range 12 < T < 16 K, whereas we 
have performed calculations at 10 K. Comparison between 
model 2 and observations toward LDN 1689N shows an 
agreement within a factor of 2 for the reported species ex- 
cept for ND2H which displays a spectacular enrichment. 
Another source of deuteration than the one assumed may 
be at work. A plausible mechanism is evaporation from 
dust grains, resulting fr om the shocks pr oduced by the 
outflows, as suggested bv lLis et alJ l 2002ah . This hypoth- 
esis is also suggested by the detection of doubly deuter- 



The results displayed for LDN 134N deserve special 
comments as the search of ND3 has been very deep and a 
very low upper limit has been obtained. The density and 
depletio n conditions are probably between th ose taken for 
model 2 llTine et all20o"fllRoueff et all200(A and model 3 
l|Pagani et al .120051) . Comparison between models and ob- 
served deuterium fractionation ratios is satisfactory ex- 
cept for ND3 where predictions are above the observed 
upper limit, and for N2H + which exhibits an impressive 
deuterium fractionation ratio towards this line of sight. 
Such values may be obtained within the present treat- 
ment if the dissociative recombination rate coefficient of 
the deuterated N2D + were significantly smaller than that 
of N2H + , a possibility which is not excluded (such effects 
have been measured for H^" and D^" but not studied ex- 
perimentally for these nitrogen bearing ions). 

The observations in the Barnard 1 molecular cloud are 
reproduced within a factor of 2 by Model 2 except for 
ND2H which is underpredicted by a factor of 3. As temper- 
atures in prestellar clouds can vary between values as low 
as 7 K in the center and ~50 K at large radii, we display 
in Figures and El the relative abundances of deuter- 
ated isotopologues of ammonia and of other molecules 
at temperatures between 5 and 50 K for the conditions 
pertaining to model 2, with a molecular hydrogen den- 
sity of 10 5 cm" 3 . The deuterium fractionation of the var- 
ious molecules is of the order of 10% for temperatures 
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between 5 and 20 K (except for ND/NH where the value 
is about unity at low temperatures), and decreases quickly 
for higher temperatures. 

The peak in the deuterium fractionation of ND3 rela- 
tive to NH 3 at T sa 12 K is directly linked to the barrier 
of the reaction of N + with H P which is of the o r der o f 
16 K l|Marauette et alJll98S|) . Indeed, IJaca et alJ l|l990|) 
measured lower ammonia fractionation ratio in hot cores 
< 0.01 than in the dense cold cores s tudied here. The 
tempe rature dependence is discussed bv lShah k, Wootter] 
ll200lh . The ammonia fractionation remains close to 10% 
for a kinetic temperature below 30 K, but drops to a few 
percent in warmer regions. The turnover in the model pre- 
sented here occurs at a slightly lower temperature (22 K), 
but the overall shape is fairly well reproduced, given the 
uncertainty in the observed temperature measurements. 
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0.001 
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Fig. 11. Results of model calculations of the ratio of dif- 
ferent deuterated isotopologues of ammonia as a function 
of temperature. The conditions of model 2 have been em- 
ployed. The full line shows NH2D relative to NH3 the 
dotted line shows ND2H relative to NH2D and the dash- 
dotted line shows ND 3 reltive to ND 2 H . Note that ND 2 H 
relative to NH2D is the lowest ratio for all the tempera- 
tures considered in this calculation. 

Another interesting feature is seen in Figure IT2l where 
the modelled DCN/HCN ratio is found to be smaller than 
that of DNC/HNC at temperatures below 20 K, whereas 
the opposite is found at higher temperatures. This trend 
is clearly found towards LDN 134N which is known to be 
a cold s ource. Such a tendency had already been empha- 
sized by iRoueff fc Gerinl l(2003h for a sample of dark cold 
sources. We have been led to carefully consider the vari- 
ous channels leading to deuterated HCN and HNC which 
are produced, in such gas phase chemistry, via the disso- 
ciative recombination of HCND+, DCNH+ and HDNC+. 
Table El displays the values assumed. 
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0.001 




Fig. 12. Results of model calculations of the ratio of 
different deuterated isotopologues of nitrogen contain- 
ing molecules as a function of temperature. The condi- 
tions of model 2 have been employed. The full line dis- 
plays ND/NH , crosses and diamonds on the dotted line 
shows DCN/HCN and DNC/ HNC respectively whereas 
the dashed-dotted line without and with open circles refer 
to DCO+/HCO+ and N 2 D+ and N 2 H+ respectively. 

Table 12. Assumed reaction rate coefficients and branch- 
ing ratios in dissociative recombination reactions of 
deuterated HCNH+ and H 2 NC+ 



Reaction 



k (in 10" 



_3„-1n 



HCND+ 


+ electron 


-» HCN + D 


1.16 (T/300)" 


-U.5 






DNC + H 


2.33 (T/300)" 


-0.5 


DCNH 4 


+ electron 


-» DCN + H 


2.33 (T/300)" 


-0.5 






HNC + D 


1.16 (T/300)" 


-0.5 


HDNC 4 


+ electron 


-» HNC + D 


0.58 (T/300)" 


-0.5 






DNC + H 


1.16 (T/300)" 


-0.5 






NHD + C 


1.75 (T/300)" 


-0.5 



We note that this prediction i s not derived in the ac- 
cretion models of iRoberts et alJ l|200A . DCN/HCN and 
DNC/HNC ratios become closer in LDN 1689N where 
the temperature m ay b e around 16 K a s discussed by 
Miz uno et al. hl99Ch and lStark et alJ (l2004) . m nice agree- 
ment with the model predictions. In addition, we display 
the temperature dependence of the deuterium fractiona- 
tion ratios of the HCO + and N 2 H + molecular ions, which 
are very similar. 

The ND 3 to NH 3 ratio is highest for LDN 1544. 
However, the ammonia column density in this source may 
be underestimated because of the large beam size of the 
observations. High fractionations are also obtained for 
Barnard 1 and LDN 1689N, for which we could benefit 
from higher quality ammonia data. We therefore believe 
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that the correction for the ammonia column density is 
at most a factor of three and the ND3 to NH3 ratio in 
LDN 1544 is close to the values obtained in Barnard 1 and 
LDN 1689N. This line of sigh t is well studied a nd known 
to be a highly depleted core IjCaselli et al . 2003). The ob- 
serv ations and the results of model 3 agree within a factor 
of 3. iRoberts et alJ l|2004|) have made detailed modelling 
of LDN 1544 by introducing the density and temperature 
structure of this prestellar core in a model including ac- 
cretion processes on dust grains and an extensive chem- 
ical network. The overall predicted deuterium fractiona- 
tion is larger by about one order of magnitude than the 
present observations. However, a nice agreement between 
models and obser vations is obta i ned fo r doubly deuterated 
formaldehyde bv IRoberts et al.N2004h . 

Extremely young protos tars have been foun d in the 
Barnard 1 molecular cloud ( Hirano et al.lll99^) and the 
depletion factor is about 3 compared to TMC-1 (cf. 
iLis et alJl2002bh . Recent observations have revealed the 
presen ce of HPS. DoS fVastel et alJ fcoofl . HDCS and 
D 2 CS llMarcelino et all2005l) in this dark cloud, which de- 
serves a specific chemical analysis. The level of deuteration 
is quite high and probably involves a source of deuteration 
in addition to that provided by gas phase chemistry alone. 

5. Conclusions 

The present study reports detection of fully deuterated 
ammonia towards 5 lines of sight, with derived ND3/NH3 
ratios between 2 xlO -4 and 1-2 xlO -3 in various pre- 
stellar cores. A very low upper limit for this ratio (1.6 
xl0~ 4 ) has been found towards the dense dark cloud 
LDN 134N. These observations imply a tremendous en- 
hancement factor over the cosmic deuterium/hydrogen 
value which varies from about 11 to over 12 orders of mag- 
nitude. The spatial extension is limited and the ND3 peak 
does not coincide with the DCO + peak in LDN 1689N, a 
source where such a study is possible thanks to its excep- 
tionally high deuterium enhancement. Deuterated species 
do not peak towards protostars themselves, but at offset 
positions, suggesting that protostellar activity decreases 
deuteration built in the pre-stellar phase. We have studied 
the origin of the deuterium enrichment in terms of a steady 
state gas phase chemical model, which includes a limited 
number of molecules, but which allows full deuteration of 
basic important intermediates. The role of grains is sim- 
ulated by introducing different values of the depletion of 
C and O so that the important deuterated molecular ions 
H 2 D+, D 2 H+, D 3 +, CH 2 D+, CD 2 H+, CD 3 + are not de- 
stroyed by reactions with O, CO, ... in the highly depleted 
regions. Our models make the following predictions: 

— Higher density models with higher C and O depletions 
yield greater enhancement of the relative abundances 
of the deuterated ammonia isotopologues. 

— The models predict abundance ratios of the deuterated 
ammonia isotopologues, which agree reasonably well 
with new and existing observations. 



The fractionation ratios of ammonia deuterated iso- 
topologues remain large for temperatures as high as 
20 K. 

The temperature dependence of the ratios DCN/HCN 
and DNC/HNC is not identical. Whereas DNC/HNC 
is larger at low temperatures, the opposite becomes 
true at higher temperatures. 

A very high deuterium fractionation is predicted for 
ND. We also find a significant abundance of NHD in 
our models. Both species have rotational transitions 
in the submillimctcr domain and are potentially de- 
tectable from the ground. 

D 2 may be an important component of the cloud and 
the corresponding reactions have to be introduced. 
There is, unfortunately, no hope of detecting it under 
the conditions which characterize cold, dark clouds. 
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